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Abstract 

' Using a new method and additional (conditional and partial) equivalence transforma- 

. tions, we performed group classification in a class of variable coefficient (1 + l)-dimensional 

' nonlinear diffusion-convection equations of the general form f{x)ut = {D{u)ux)x + K{u)Ux- 

. We obtain new interesting cases of such equations with the density / localized in space, which 

' have non-trivial invariance algebra. Exact solutions of these equations are constructed. We 

I ■ also consider the problem of investigation of the possible local transformations for an arbi- 

trary pair of equations from the class under consideration, i.e. of describing all the possible 
partial equivalence transformations in this class. 
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The problems of group classification and exhaustive solutions of such problems are interesting 



^ ■ 1 Introduction 
in 
cn 

\^ . not only from purely mathematical point of view, but is also important for applications. In 

\ physical models there often exist a priori requirements to symmetry groups that follow from 

■ physical laws (in particular, from Galilei or relativistic theory). Moreover, modelling differential 

. equations could contain parameters or functions which have been found experimentally and so 

p, ^1 are not strictly fixed. (It is said that these parameters and functions are arbitrary elements.) At 

the same time mathematical models have to be simple enough to analyze and solve them. Solving 
the problems of group classification makes possible to accept for the criterion of applicability 
the following statement. Modelling differential equations have to admit a group with certain 
properties or the most extensive symmetry group from the possible ones. 

In this paper we consider a class of variable coefficient nonlinear diffusion-convection equations 
^ , of the form 

. . f{x)ut = {g{x)D{u)u,)^ + K{u)u^, (1) 

where / = /(x), g = g{x), D = D{u) and K = K{u) are arbitrary smooth functions of their 
variables, f{x)g{x)D{u) 7^0. The linear case of (1) {D,K = const) was studied by S. Lie [1] in 
his classification of linear second-order PDEs with two independent variables. (See also modern 
treatment of this subject in [2].) That is why we assume below that {Du,Ku) 7^ (0,0), i.e. (1) 
is a nonlinear equation. 

Moreover, using the transformation i = t, x = J u = u, we can reduce equation (1) to 

/(i)% = {D{u)ui)s: + K{u)ui, 

where f{x) = g{x)f{x) and g{x) = 1. (Likewise any equation of form (1) can be reduced to 
the same form with f{x) = 1.) That is why without loss of generality we restrict ourselves to 
investigation of the equation 

f{x)ut = iD{u)u,)^ + Kiu)u^. (2) 
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Apart from their own theoretical interest, equations (2) are used to model a wide variety of 
phenomena in physics, chemistry, mathematical biology etc. For the case /(x) = 1 equation (2) 
describes vertical one-dimensional transport of water in homogeneous non-dcformable porous 
media. When K{u) = this equation describes stationary motion of a boundary layer of fluid 
over a flat plate, an vortex of incompressible fluid in porous medium for polytropic relations 
of gas density and pressure. The outstanding representative of the class of equations (2) is the 
Burgers equation that is the mathematical model of a large number of physical phenomena. (For 
more detail refer to [4-7].) 

Investigation of the nonlinear heat equations by means of symmetry methods started in 1959 
with the Ovsiannikov's work [8] where he studied symmetries of the equation 

ut = {f{u)ux)x- (3) 
In 1987 I.Sh. Akhatov, R.K. Gazizov and N.Kh. Ibragimov [9] classified the equations 

Ut = G{ux)uxx- (4) 
V.A. Dorodnitsyn (1982, [10]) performed group classification of the equation 

Ut = {G{u)ux)x + 9{u). (5) 

A. Oron, P. Rosenau (1986, [11]) and M.P. Edwards (1994, [12]) presented the most extensive 
at the time list of symmetries of the equations 

Ut = {G{u)u.j,)-j: + f{u)Ux. (6) 

The results of [8, 10, 11] were generahzed by R.M. Cherniha and M.I. Serov (1998, [13]) who 
classified the nonlinear heat equation with convection term 

Ut = {G{u)Ua;)x + f{u)U:c + g{u). (7) 

It should be noted that equations (l)-(7) are particular cases of the more general class of 
equations 

Ut = F{t, X, U, Ux)Uxx + G{t, X, u, Ux). (8) 

Group classification of (8) is presented in [17-21]. However, since the equivalence group of (8) 
is essentially wider than those for (1)^(7) the results of [17-21] cannot be directly used to sym- 
metry classification of equations (l)-(7). Nevertheless, these results are useful to find additional 
equivalence transformations in the above classes. 

Equations of form (2) have been also investigated from other than classic Lie symmetry points 
of view. For instance, potential symmetries of subclasses of (2) where e.g. either / = 1 or = 
were studied by C. Sophocleous [14-16]. 

Some symmetry properties of class (1) were considered in a recent paper [3]. However, it does 
not present correct and complete results on the subject. A overwhelming majority of cases with 
nontrivial Lie symmetry were omitted, and there are mistakes in the cases that were adduced. 
Nevertheless the subject seemed very interesting, and we decided to study Lie symmetries of 
class (2). 

The ultimate goal of this paper is to present an example of exhaustive solution of the group 
classification problem in quite a difficult case. After giving a precise definition and a discussion 
of this problem in general case, we performed the complete extended group classification, found 
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additional equivalence transformations and exact solutions of equations (2). A lot of new inter- 
esting cases of extensions of maximal Lie symmetry group were obtained for these equations. 
For example, we determined the equations which have the density / localized in the space of 
X and are invariant with respect to m-dimensional (2 < m < 4) Lie symmetry algebras, that 
allows construction of new exact non-stationary solutions for them. 

Problems of general group classification, except for really trivial cases, arc very difficult. It 
can be illustrated by the multitude of papers where such general classification problem is solved 
incorrectly or incompletely. There are also many papers on "preliminary group classification" 
where authors list some cases with new symmetry but do not claim that the general classification 
problem is solved completely. For this reason finding an effective approach to simplification is 
essentially equivalent to feasibility of solving the problem at all. In this paper we develop and 
apply a simple and effective tool based on investigation of specific compatibility of classifying 
conditions. First it was proposed in [23] and then applied to solving a number of different group 
classification problems [24-26]. The other tool is systematically to use additional (conditional 
and partial) equivalence transformations that allows us to put in order, verify and analyze the 
obtained results. 

Our paper is organized as follows. First of all (Section 2) we describe the group classification 
method used here and introduce the notions of conditional and partial equivalence. Then (Sec- 
tions 3) we significantly enhance the results of [3] and give the complete group classification of 
class (2). Since the case f{x) = l has a great variety of applications and was investigated earlier 
by a number of authors, we collect results for this class together in Section 4. Section 5 contains 
the proof of the main theorem on group classification of the class (2). We attempted to present 
our calculations in most reasonable detail so as it would be feasible to verify that. Conditional 
equivalence transformations are considered in Section 6, where we also present four lemmas on 
possible local equivalence transformations between two arbitrary equations of form (2) . The re- 
sults of group classification are used to find exact solutions of equations from class (2) (Section 7). 



2 Group classification method and additional equivalence 

Let us describe the classical algorithm of group classification restricting, for simplicity, to the 
case of one differential equation of the form 

L\x,u^n)) = L(a;,U(„),^(a;,U(„))) = 0. (9) 

Here x = (xi, . . . ,xi) denotes independent variables, u is a dependent variable, is the set 
of all the partial derivatives of the function u with respect to x of order no greater than n, 
including u as the derivative of zero order. L is a fixed function of x, and 9. 9 denotes the 
set of of arbitrary (parametric) functions 0(a;,U(„)) = (^''^(a;, U(„)), . . . ,9^{x,U(^n))) satisfying the 
conditions 

5'(a;,U(„),6'(g)(a;,U(„))) = 0, S = {Si, . . . ,Sr). (10) 

These conditions consist of r differential equations on 9, where x and 'fX(„) play the role of 
independent variables and 0(g) stands for the set of all the partial derivatives of 9 of order no 
greater than q. In what follows we call the functions 9{x,u^n)) arbitrary elements. Denote the 
class of equations of form (9) with the arbitrary elements 9 satisfying the constraint (10) as L\s. 

Let the functions 9 be fixed. Each one-parametric group of local point transformations that 
leaves equation (9) invariant corresponds to an infinitesimal symmetry operator of the form 

Q = ^"(x, u)dx^ + r]{x, u)du. 
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(here the summation over the repeated indices is understood) . The complete set of such groups 
generates the principal group G^^^ = G^^^{L,9) of equation (9). The principal group G™*^^ 
has a corresponding Lie algebra A^^^ = A^^^^^{L,9) of infinitesimal symmetry operators of 
equation (9). The kernel of principal groups is the group 

0:S{e)=O 

for which the Lie algebra is 

e:S{d)=0 

Let G^'i^" = G^"^^" {L, S) denote the local transformations group preserving the form of equa- 
tions from L\s. (Sometimes one consider a subgroup instead the complete equivalence group.) 

The problem of group classification consists in finding of all possible inequivalent cases of 
extensions of A"^'^^, i.e. in a listing all G ^^"'^-inequivalent values of 9 that satisfy equation (10) 
and the condition A'°^^(L,6') 7^ A^"'. 

In the approach used here group classification is application of the following algorithm [2, 22] : 

1. From the infinitesimal Lie invar iance criterion we find the system of determining equations 
for the coefficients of Q. It is possible that some of the determining equations does not 
contain arbitrary elements and therefore can be integrated immediately. Others (i.e. the 
equations containing arbitrary elements explicitly) are called classifying equations. The 
main difficulty of group classification is the need to solve classifying equations with respect 
to the coefficients of the operator Q and arbitrary elements simultaneously. 

2. The next step involves finding of the kernel algebra A^^^ of principal groups of equations 
from L\s. After decomposing the determining equations with respect to all the uncon- 
strained derivatives of arbitrary elements one obtains a system of partial differential equa- 
tions for coefficients of the infinitesimal operator Q only. Solving of this system yields the 
algebra A^^^. 

3. In order to construct the equivalence group G^"^^^^ of the class LI5 we have to investigate 
the local symmetry transformations of system (9), (10), considering it as a system of partial 
differential equations with respect to 9 with the independent variables x, tt(„). Usually 
one considers only transformations being projectible on the space of the variables x and 
u. Although in the case when 9 depends, at most, on these variables, it can be assumed 
the transformations of them depend on 9 too. After restricting ourselves to studying of 
the connected component of unity in G^'^^^'^, we can use the Lie infinitesimal method. To 
find the complete equivalence group (including discrete transformations) we are supposed 
to use a more complicated direct method. 

4. If A^^ is an extension of ^4^'''' (i.e. when ^™^^(L, 9) ^ A^^^), then the classifying equations 
define a system of nontrivial equations for 9. Depending on their form and number, we 
obtain different cases of extensions of ^4'^'^''. To completely integrate the determining equa- 
tions we have to investigate a large number of such cases. In order to avoid cumbersome 
enumeration of possibilities in solving the determining equations we can use, for instance, 
algebraic methods [17-21,27], a method which involves the investigation of compatibility 
of the classifying equations [23-26] or combined methods [28,29]. 
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The result of application of the above algorithm is a list of equations with their Lie invariance 
algebras. The problem of group classification is assumed to be completely solved if 

i) the list contains all the possible inequivalent cases of extensions; 

ii) all the equations from the list are mutually inequivalent with respect to the transformations 

from G^^^'^; 

in) the obtained algebras are the maximal invariance algebras of the respective equations. 

Such list may include equations being mutually equivalent with respect to local transformations 
which do not belong to G^'i™^. Knowing such additional equivalences allows to simplify essen- 
tially further investigation of L\s. Constructing them can be considered as the fifth step of the 
algorithm of group classification. Then, the above enumeration of requirements to the resulting 
list of classification can be completed by the following step: 

iv) all the possible additional equivalences between the listed equations are constructed in 
explicit form. 

One of the ways for finding additional equivalences is based on the fact that equivalent equa- 
tions have equivalent maximal invariance algebras. The second way is systematical study of 
conditional and partial equivalence transformations in the class L\s. Let us give definition of 
such transformations. Consider a system 



G^'^^"{L, {S,S')) denote the equivalence group of the subclass L\s,s' of L\s, where the func- 
tions 9 satisfy systems (10) and (11) simultaneously. 

Notion 1. We call the transformations from G^'^^^'^{L, {S, S')) a conditional equivalence trans- 
formations of class L\s (under the additional constraint S' ). The local transformations which 
transform equations from L\s^s' io L\s are called partial equivalence transformations of the 
class L\s (under the additional constraint S' ). 

It is obvious that any conditional equivalence is a partial one under the same additional 
constraint and any local symmetry transformation of equation (9) for a fixed value 9 = 9^{x, ■«(„)) 
is a partial equivalence transformation under the constraint 9 = 9^. The problem of description 
of all the possible partial equivalence transformations in the class L\s is equivalent to that 
of local transformations between two arbitrary equations from L\s. Additional constraints on 
arbitrary elements never imply constraining sets of partial equivalences. 

3 Results of classification 

Consider a one-parameter Lie group of local transformations in {t,x,u) with an infinitesimal 
operator of the form Q = ^*(t, x,ii)9i + S,^{t,x,u)dx + r]{t, x,u)du, which keeps equation (2) 
invariant. The Lie criterion of infinitesimal invariance yields the following determining equations 
for and rj 




= f * = f ^ = 

Si Su Su 

Driuu + D^rju - D^{2C - it) + DuuV - y = 0, 
2C - it + y r = ^V, fVt - Kr,, - D^,, = 0, 
K (fye + C - it) + mix - 2%«) - '^DuVx - KuV - f^t = 0- 




(12) 
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Investigating the compatibility of system (12) we obtain an additional equation rjuu = without 
arbitrary elements. With this condition system (12) can be rewritten in the form 

C = ^i = C = 0, Vuu = 0, (13) 

2€ - et + ^-fe = ^v, (14) 

Dvxx + Kr], - fru = 0, (15) 

(D^K - K^D) ^-KC- 2D^Vx + D^l, - /C? - 2D7?,„ = 0. (16) 
Equations (13) do not contain arbitrary elements. Integration of them yields 



Thus, group classification of (2) reduces to solving of classifying conditions (14)-(16). 

Splitting of system (14)-(16) with respect to the arbitrary elements and their non-vanishing 
derivatives gives the equations = 0, = 0, ?7 = on the coefficients of operators from A^^^ 
of (2). As a result, the following theorem is true. 

Theorem 1. The Lie algebra of the kernel of principal groups of (2) is A^'^'^ = {dt). 

The next step of algorithm of group classification is finding of equivalence transformations 
of class (2). To find these transformations, we have to investigate Lie symmetries of the system 
that consists of equation (2) and additional conditions 

ft = fu = 0, Dt = D^ = 0, Kt = K^ = 0. 

Using the classical Lie approach, we find the invariance algebra of the above system that 
forms the Lie algebra of G*^^"'"^ for class (2). Thus, we obtain the following statement. 

Theorem 2. The Lie algebra of G'^'^''' for class (2) is 

^equiv ^ l^Q^^ Q^^ Q^^ ^ fQ^^ _ 2fdf - KdK, Udu, fdf + KdK + DOd). (18) 

Therefore, G^^^^'^ contains the following continuous transformations: 
i = te^* + £i, X = xe^^ + £2, u = ue^^ + £3, 

where £i,...,e7 are arbitrary constants. For class (2) there also exists a nontrivial group of 
discrete equivalence transformations generated by four involutive transformations of alternate 
sign in the sets {t, D, K}, {x, K}, {u} and {/, D, K}. It can be proved by the direct method that 
Qequiv coincides with the group generated by the both continuous and discrete transformations 
from the above list. 

Theorem 3. A complete set of inequivalent equations (2) with respect to the transformations 
from G'^^^" with A^^ ^ A^*^^ is exhausted by cases given in Tables 1-3. 
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Table 1: Case \fD{u) 



N 


K{u) 


fix) 


T) • J? (1 max 

Basis 01 A 


1 


V 


V 


dt 


2a 


V 


sx 

e 


dt, etdt + 


2b 


L> 




ft, -ftdt - e^'dx 


2c 


D 




dt, + 2)tdt + '^)e^dx 


2d 







dt, {v + 2)tdt + xdx 


2e 


1 




dt, e~*{dt — xdx) 


3a 





1 


dt, dx, 2tdt + xdx 


3b 


D 




dt, 2tdt-dx, e'dx 



Here7,i/f^0, e = 0, 1 mod C^""'^, 7 = ±1 mod G^""'^. 
Additional equivalence transformations: 

1. 2b ^ 2a (if = 0, e = 1): t = t,x = 76"'", u = u; 

2. 2c (v / -2) -^2&{K = -D/{v + 2), e = 1): f = t, 5 = (i^ + 2) In je"^ + 7|, m = u; 
2c {v = -2) -^2&(K = -D, e = 0): i=t,x = In \e~'' +^\,u = u; 

3. 2d {f ^ -2) ^ 2a (is: = --D/(i' + 2), £ = 1): i = t, x = + 2)ln\x\, u = u; 
2d (z/ = -2) ^ 2a (if = -D, e = 0): f = t, x = In |a;|, u = u; 

4. 2e ^ 2a (is: = -D, e = 1): t = e*, x = In |x| + 1, m = w; 

5. 3b — » 3a: i = t, x = e~^, m = w. 



Table 2: Case D{u) = e^^" 



N 




is:(u) 


fix) 


Basis of yl""^" 


1 


V 




\x\^ 


ft, (A^ — Aj/ + — 2!y)tft + (p - + 9u 


2 


V 


e" 


1 


ft, dx, in - 2)tdt + (m - + du 


3 


1 


u 


1 


dt, dx, tdt + ix- t)dx + du 


4 


1 


ee" 


V 


dt, tdt - du 


5a 


1 





f\x) 


dt, tdt — du, atdt + ilSx"^ + 71a: + 7o)ft + jSxdu 


5b 


1 


u 

e 


fix) 


dt, tdt - du, atdt - iPe-'' + 71 + 7oe^)a^ + Pe'^'du 


5c 


1 


1 




dt, xdx + du, e~\dt - xdx) 


6a 


1 





1 


dt, tdt - du, 2tdt + xdx, dx 


6b 


1 


e" 




dt, tdt-du, 2tdt-dx, e^dx 


6c 


1 


e" 


x-^ 


dt, tdt-du, (e-" + 7)e"ft + 9„, -(e"" + 7) Va. + (e"" + 7)9„ 


6d 


1 





ft, tdt — du, xdx — du, x^dx + xdu 



Here A 7^ 0, £ € {0, 1} mod G"'^"'"', a, (3, 71, 70 = const and 

\ f —SPx — 271 + a 



f {x) = exp 



■ dx 



f {x) = exp 



/3e-" - 27oe" 



/9a;2 + 71a; + 70 

Additional equivalence transformations: 

1. 5b ^ 5a: i = t, x = e~^, u = u; 

2. 5c ^ 5a (a = 70 = 1, /? = 71 = 0, = x~^): i = e*, x = e*x, u = u 

3. 6b — » 6a: i = t, x = e~^, u = u; 

4. 6c ^ 6a: i= tsign(e~'" + 7), 5 = l/(e"'' + 7), m = m - In je""^ + 7]; 

5. 6d — » 6a: i = t sign a;, x = 1/x, « = w — In |a;|. 



■ dx 
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Table 3: Case D{u) = u^" 



N 




K{u) 


fix) 


Basis of 


1 


V 


u" 




5t, (/i + A/i - 2i/ - A!/)t9t + (yu - i^)a;9a, + 


2 


V 


u" 


1 


9t, dx, ip — 2v)tdt + ip~ y)xdx + udu 


3 


V 


Inti 


1 


dt, dx, Htdt + inx -t)dx + udu 


4 


V 


ew" 


V 


dt, ntdt — udu 


5a 


V 





fix) 


dt , fJ-tdt — udu , 

atdt + ((1 + p)Px^ + 71 a; + -yo)dx + 0xudu 


ob 


V 




J ix) 


dt, jJLtdt — udu, 

atdt - ((1 + M)/3e-== + 71 + 7oe==)aa, + pe-^'udu 








Jj 


Pi ^ d ^ t Pt , 'y^Pt \ it'yf~t 1 ^ 1 Pt 

ut, c yyt ^fJx}, fJ^UyUx ^ UjUu 


6a 


M -4/3 





1 


dt, ntdt — udu, dx, 2tdt + xdx 


6b 


/i / -4/3 






dt, fj.tdt — udu, "itdt - dx, e'dx 


6c 


-1 







dt, tdt + udu, dx — 'fudu, 2tdt + xdx - ^xudu 


6d 


-1 






dt, tdt + udu, e-^dx + ludu, 2tdt - dx - 'ye~'^udu 


6c 


/i ^ -4/3, -1 





\x\ 1+f 


dt, ntdt - udu, (2 + ij)tdt - (1 + n)xdx, 

(1 + fj.)x'^dx + xudu 


6f 


/i 7^ -4/3, -1 


u'^ 




dt, fitdt - udu, (2 + fi)tdt + (1 + M)(e"" + 7)e"9x, 
-(1 + M)(e-" + -yfe'^dx + (e"" + -f)udu 


(e-«'+7) l+A" 


6g 


-3/2 


1 




dt, e~*{dt — xdx), Sxdx — 2udu, e*{x^dx — 2xudu) 


7a 


-4/3 





1 


dt, 4:tdt + iudu, dx, 2tdt + xdx, x^dx - Sxudu 


7b 


-4/3 






dt, Udt + 3udu, 2tdt - dx, e-^{dx + iudu), e^dx 


8 





u 


1 


dt, dx, 2tdt + xdx - udu, tdx - du, 
t^dt + txdx — itu + x)du 



Here /x for Cases 4-6. e = 0, 1 mod G "'^'^^^ , A ^ 0, a, /3, 71 , 70 = const and 



/ (x) = exp 



-(4 + 3/i)/3x-27i + a 



dx 



f (x) = exp 



(1 + fi)l3x'^ + 71a; + 70 

Additional equivalence transformations: 

1. 5b 5a: t = t, x — e~^, u — u; 

2. 5c — > 5a (a = 70 = 1, /3 = 71 = 0, = x~^): i = e*, x = e'x, u = u 

3. 6b — » 6a: t = t,x = e~^, u = u; 

4. 6c ^ 6a (/i = —1): t = t, x = x, u = e^^u; 

5. 6d — »■ 6a (/X = —1): t = t, x = e~^, u = e'" u; 

_ 1 

6. 6e ^ 6a: t = t, x = —1/x, u = \x\ i+z'u; 

7. 6f ^ 6a: i=t,x = -l/ie'" + 7), w = le"^ + 7|"Ww; 

8. 6g ^ 6a (/X = -3/2): i=e*,x = -e~*/x, u = \e*x\~^u; 

9. 7b — »• 7a: i = t, x = e"'^, u = u. 



(2 + At)/?e-^ - 27oe"^ - a 
(1 + ^i)/3e-^ + 71 +7oe'» 



dx 
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In Tables 1-3 we list all possible G^'i"''^-inequivalent sets of functions f{x), D{u), K{u) and 
corresponding invariancc algebras. Numbers with the same arable figures correspond to cases 
that are equivalent with respect to a local equivalence transformation. Explicit formulas for 
these transformations are adduced after the Tables. The cases numbered with different arable 
figures are inequivalent with respect to local equivalence transformations. In order to simplify 
the presented results, in the case f{x) = 1 we just use the conditional equivalence transformation 
X = X — et, K = K + e (the other variables are not transformed) from (see Section 4). 

Other conditional equivalence transformations are considered in Section 6. 

Below for convenience we use double numeration T.N of classification cases and local equiv- 
alence transformations, where T denotes the number of table and N does the number of case 
(or transformation) in Table T. The notion "equation T.N" is used for the equation of form (2) 
where the parameter-functions take values from the corresponding case. 

The operators from Tables 1-3 form bases of the maximal invariance algebras iff the cor- 
responding sets of the functions f,D,K are G '^'I'^'^-inequivalent to ones with most extensive 
invariance algebras. For example, in Case 3.1 (/U, i^) 7^ (0,0) and A 7^ —1 if = 0. And in 
Case 3.2 (/x, v) ^ {(—2, —2), (0, 1)} and v ^ 0. Similarly, in Case 2.1 the constraint set on the 
parameters u and A coincides with the one for Case 3.1, and we can assume that = 1 or 
ly = 1. In Case 2.2 we consider v = 1 immediately. 

After analyzing the obtained results, we can state the following theorem. 

Theorem 4. If an equation of form (2) is invariant with respect to a Lie algebra of dimension 
no less than 4 then it can be reduced by means of local transformations to one with f{x) = 1. 



4 Group classification for subclass with f{x) = 1 

Class (2) includes a subclass of equations of the general form 

ut = {D{u)u^)^ + K{u)u^. (19) 

(i.e. the function / is assumed to be equal to 1 identically). Symmetry properties of equa- 
tions (19) were studied in [11, 12]. But we are not aware of any paper containing correct and 
exhaustive investigation on the subject. Now let us single out the results of group classification 
of equations (19) from the above section. 

Theorem 5. The Lie algebra of the kernel of principal groups of (19) is = {dt,dx)- 
Theorem 6. The Lie algebra of the equivalence group Q^"^^" for the class (19) is 

Ai^q^i^ = {dt, dx, du, udu, tdx - Ok, 2tdt + xdx - KdK, tdt - Ddo - KOk). (20) 
Any transformation from Q'^'^^" l^as the form: 



i = teles + El, X = x£4 + ert + e2, u = uee + £3, 
D = \ k = Kel^e-^^ - £7, 

where ei, . . . , £7 are arbitrary constants, e^e^e^ 7^ 0. 



(21) 



Theorem 7. The complete set of G^^^^^ -inequivalent extensions of for equations (19) is 
given in Table 4- 
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Table 4: Case /(x) = 1. 



N 


Din) 




Basis of A"""" 


1 


V 


V 


ft, 




2 


V 





dt, dx, 


2tdt + xdx 


3 




e" 


dt, dx, 


in - 2)tdt + (a* - l)xdx + du 


4 


e" 


u 


dt, dx, 


tdt + (a: — t)dx + du 


5 


e" 





dt, dx, 


tdt - du, 2tdt + xdx 


6 


u'^ 




dt, dx, 


{ji — 2v)tdt + (m - v)xdx + udu 


7a 







dt, dx, 


fitdt - udu, 2tdt + xdx 


7b 






dt, dx, 


2tdt + udu,er'^'{dx + udu) 


8 







dt, dx, 


Atdt + 3udu, 2tdt + xdx, x^dx — ixudu 


9 




Inw 


dt, dx, 


IMtdt + (/ix — t)dx + udu 


10 


1 


u 


dt, dx, 


t^dt + txdx — {tu + x)du, 2tdt + xdx — udu, tdx — du 



Here ij,,^ = const, (/i, v) =^ (—2, —2), (0, 1) and for N = 6. —4/3 for N = 7a. Case 7b can be reduced 
to 7a (/u = —2) by means of the conditional equivalence transformation i =t, x = , u = e~^u. 



5 Proof of Theorem 3 

Our method is based on the fact that the substitution of the coefficients of any operator from 
^max\^^ker -j^^^ classifying equations results in nonidentity equations for arbitrary elements. 
In the problem under consideration, the procedure of looking for the possible cases mostly 
depends on equation (14). For any operator Q G A"^^^ equation (14) gives some equations on D 
of the general form 

{au + b)Du = cD, (22) 

where a,b,c = const. In general for all operators from A"^^^ the number k of such independent 
equations is no greater then 2 otherwise they form an incompatible system on D. is an 
invariant value for the transformations from G^qu'v Therefore, there exist three inequivalent 
cases for the value of k: k = 0, k = 1, k = 2. Let us consider these possibilities in more detail, 
omitting cumbersome calculations. 

fc = (Table 1). Then the coefficients of any operator from A™^ are to satisfy the system 

^ = 0, 2C-et + je = 0, -KC + Dil,-J^t = Q. (23) 

Let us suppose that if ^ (1, D).l\, follows from the last equation of the system (23) that 
Q — — 0- Therefore, the second equation is a nonidentity equation for / of the form = nf 
without fail. Solving this equation yields Case 2a. 

Now let K G (1, D), i.e. K = eD + /3 where e G {0, 1}, P = const . Then the last equation 
of (23) can be decomposed into the following ones 

The equation {S,^{fx/ f + ^e))^ = is a differential consequence of the reduced determining 
equations. Therefore, the condition fx/ f + 2e = is a classifying one. 

Suppose this condition is true, i.e. / = e~^^^ mod G^'^^^'^. There exist three different possi- 
bilities for values of the parameters e and 

£ = 1, e = l, f3 = 0; £ = (then /3 = mod G^'^''"), 
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which yield Cases 2a, 3b and 3a respectively. 

Let e = and fxjj 7^ 0. Then either our consideration is reduced to Case 2a or / = 
mod G^'^^" where i-i ^ 0. Depending on the value of the parameter /? (/3 = or /? 7^ and then 
/J, = —1) we obtain Case 2d or Case 2e. 

Let £ = 1 and fx/f 7^ —2. Then (3 = and fx/f = (Cie^ + Co)~^ — 2 where we assume 
Ci 7^ to exclude Case 2a. Integrating the latter equation depends on whether Cq vanishes or 
not, and results in Cases 2b and 2c. 

fc = 1. Then D G {e", u'', /i 7^ 0} mod G^^"^^ and there exists Q G with 7^ 0. 

Let us investigate the first possibility D = (Table 2). Equation (14) implies rju = 0, i.e. 
r] = r){t,x). Therefore, equations (14)-(16) can be written as 

2C - i\ + ^-fe = V, e-rixx + Krix - frit = 0, 

/ (24) 

{K - Ku)v - m - + e^i^xx - 2%) = 0. 

The latter equation looks like = uK + 6e" + c with respect to if, where v,h,c = const. 
Therefore, K is to take one of the following five values. 

1. K = e^" + xie" + x-Q mod G'^'^^" ^ where 7^ 0, 1. (Here and below Xj = const, i = 0, 1.) 
Then rj = const, xi = 0, and either xq = if / 7^ const or xq = mod Gl'^™^ if / = const that 
imply = 0, = 0, therefore / = \x\^ mod G^^"'^ (Cases 1 and 2). 

2. K = n + xie" + xq. In the way analogous to the previous case we obtain xi = 0, / = 1 
mod G^'i"'^, xo = mod G^''''' (Case 3). 

2. K = ue^ + xie" + xq mod G^^^'^ . It follows from system (24) that ry = for any operator 
from A™^, i.e. we have a contradiction to assumption r/ 7^ for some operator from A^^. 

4. K = + XQ. Then r/^ = C{t)e-'' + (^{t), = cj^(t)e^ + cT°(t) - C{t)e-''. It can be proved 
that C,l = = = = 0, cither xq = if / 7^ const or xq = mod G^^^" if / = const, and 
therefore = 0. The first equation of (24) implies that the function / is to satisfy / (/ = 0, 1, 2) 
equations of the form 

fx /3e"^ - a- 



f /3e ^ + 71 + 7oe^ 



with non-proportional sets of constant parameters (a, /3, 70, 71). The values / = and / = 1 
correspond to Cases 4 (e = 1) and 5b. An additional extension of exists for Z = 2 in 

comparison with Z = 1 iff / is a solution of the equation 

fx ^ Aae-^ 

/ Aie-^ + Ao ' 

where either A2 = or A2 = 3Ai 7^ 0. Integrating the latter equation gives Cases 6b and 6c. 

5. K = XQ. Then r?^ = (^(t)x + C^it), = a^it)x + a'^{t) + C^{t)x'^. It follows from compatibility 
of system (24) that rjt = = ii f ^ {x'^, 1} mod G^'i"^^ or xq = 0. The values / = 
Xq 7^ result in Case 5c. If / 1} mod G^"^^^"^ and xq = 0, we obtain Case 1 with u = 0. 

If / = const then xq = mod G^"^^^^ . Below xq = 0. The first equation of (24) holds when the 
function / is a solution of a system of I (I = 0, 1, 2) equations of the form 

fx _ -3/3x + g - 271 
/ + 71X + 70 
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with non-proportional sets of constant parameters (a, /3, 70, 71). The values I = and I = 1 
correspond to Cases 4 (e = 0) and 5a. Additional extensions for Z = 2 exist iff / is a solution of 
the equation 

fx A2 



/ Aix + Ao 

where either A2 = or A2 = — 3Ai 7^ 0. These possibilities result in Cases 6a and 6d. 

Consider the case D = u^' (Tables 3). Equation (14) implies 77° = 0, i.e. r] = r]^{t,x)u. 
Therefore, system (14)-(16) can be written as 

K - + je = f^v\ n'^vL + Kvl - fvl = 0, ^25) 

{^^K - uK^W - + {ilx - 2(/x + lW.)u^^ - fCf = 0. 

The latter equation looks with respect to K similarly to uK^ = uK+hu^+c, where i',b,c = const. 
Therefore, K is to take one of the five values. 

1. K = u'^ + Kiu^ + xq mod where v 7^ 0, /U. Equations (25) imply ry^ = const, 
= {fi — u)'q^x + (T(t), xiQ = (therefore, xi = since t?^ = 0), / = mod Q^'^^^'' , 

= (^ + A/i - 2z/ - Al/)ry^ \a = 0, and either xq = if A / (Case 1) or xq = mod G^"^'^ if 
A = (Case 2). 

2. K = Inu + x\u^^ + xo mod G^^^"" . In the way analogous to the previous case we obtain 
XI = 0, / = 1 mod G^i"'^, xo = mod G^^'i"''' (Case 3). 

3. if = u^^\nu + xiu^ + HQ mod G^"^™^ . It follows from system (25) that = for any operator 
from yl™^^^ i.e. we have a contradiction to assumption r] ^ for some operator from A^^^. 

4. K = n^^ + xo mod G^""'^. Then r/^ = C^{t)e-'' + C^{t), = (Ti(t)e^ + (T°(t) - (^+ l)C^(t)e-^. 
It can be proved that (^^ = Ct* = cr^ = = 0, either xq = if / 7^ const or xq = mod G^*^"'^ 
if / = const, and therefore = 0. The first equation of (25) implies that the function / is to 
satisfy I (Z = 0, 1, 2) equations of the form 

^ ifi + 2)/3e-'^ -a- 2706^ 
/ (m + l)/3e-^ + 71 + 7oe^ 

with non-proportional sets of constant parameters (a, /?, 70, 71). The values I = and I = 1 
correspond to Cases 4 (e = 1) and 5b. Z = 2 iff / is a solution of the equation 

^ A2e-^ ^ 
/ Aie-^ + Ao 

Looking for the inequivalent possibilities of integrating of this equation results in Cases 6b, 6d, 
6f, 7b. 

5. K = XQ. Then r]^ = CHt)x + C°(i), = (^\t)x + a^{t) + (n + l)CHt)x'^- It follows from 
compatibility of system (25) that 774 = = if / ^ {x"^, 1} mod G^"'"^ or xq = 0. The values 
/ = x~^, Xo 7^ result in Cases 5c and 6g. If / ^ {x~^, 1} mod G^'^^^'^ and xq = 0, we obtain 
Case 1 with 1/ = 0. li f = const then xo = mod G^'^"^^. Below xo = 0. The first equation 
of (25) holds when the function / is a solution of a system of / (Z = 0, 1, 2) equations of the form 

^ ^ ~{3fi + 4)Px + a-2-fi 
f {11+ l)/3x2 + 7ix + 70 
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Table 5: Conditional equivalence algebras 



Conditions 


Basis of 


K = D 
K = D = e^ 
D = e^,K = Q 
D = K = u'^ 
D = u'',K = 


dt, d., du, udu, tdt + fdf, e"(a. - 2/9/), fdf + Ddo 

dt, tdt + fdf, dx, du + fdf, xdx - "2 fdf, x^dx + xdu - 3xdf 

dt, tdt + fdf, dx, du + fdf, xdx-2fdf, x'^dx + xdu - 3xdf 

dt, tdt + fdf, dx, du + nfdf, e^idx - 2/9/), e-''{{l+n)dx - udu + (2+/x)/9/) 

dt, tdt + fdf, dx, du+ufdf, xdx-2fdf, {1 + ij.)x'^dx + xudu - + Sn)xfdf 



with non-proportional sets of constant parameters (a, /3, 70, 71). The values / = and / = 1 
correspond to Cases 4 (e = 0) and 5a. Z = 2 iff / is a solution of the equation 

fx _ A2 

/ Xix + Xo' 

Looking for the inequivalent possibilities of integrating of this equation results in Cases 6a, 6c, 
6e, 7a. 

k — 2. Assumption on two independent equations of form (22) on D yields D = const, i.e. 
D = 1 mod G^'^^^'^. Ku 7^ (otherwise, equation (2) is linear). Equations (14)-(16) can be 
written as 

2^. - il + y r = 0, + Kt), - flit = 0, 

-K^n-Kil + e,,-fit-2r,l = Q. 

The latter equation looks similarly to {au + b)Ku = cK + d with respect to K, where a, b,c,d = 
const. Therefore, to within transformations from G^'^^" K is to take one from four values: 

K = + xq, V ^ K = lnu + xo; K = e" + xo; K = u. 

Classification for these values is carried out in the way similar to the above. The obtained 
extensions can be entered in either Table 2 or Table 3. 

The problem of the group classification of equation (2) is exhaustively solved. 

6 Additional equivalence transformations 

When we impose some restrictions on arbitrary elements we can find additional equivalence 
transformations named conditional equivalence transformations (see Notion 1). As mentioned 
above, the simplest way to find such equivalences between previously classified equations is 
based on the fact that equivalent equations have equivalent maximal invariance algebras. A 
more systematic way is to classify these transformations using the infinitesimal or the direct 
methods. Examples of conditional equivalence algebras calculated by the infinitesimal method 
are listed in Table 5. 

To find the complete collection of additional local equivalence transformations including both 
continuous and discrete ones, we are to use the direct method. Moreover, application of this 
method allows us to describe all the local transformations that are possible for pairs of equations 
from the class under consideration. The problem of this sort was first investigated for wave 
equations by Kingston and Sophocleous [30-32]. Now we state a number of simple but very 
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(26) 



useful lemmas containing preliminary results on solution of this problem. (We imply that the 
condition of nonsingularity is satisfied for all the transformations.) 

Lemma 1. [33] For any local transformation between two evolutionary second-order equations 
(i.e. equations of the form ut = H(t,x,u,Ux,Uxx) where H^^^ ^) ^he transformation of the 
variable t depends only on t. 

Lemma 2. Any local transformation between two evolutionary second-order quasi-linear equa- 
tions having the form ut = F{t, x, u)uxx + G{t, x, u, Ux) where F ^ is projectible, i.e. t = T(t), 
X = X(t, x), u = U{t, X, u). 

Lemma 3. Any local transformation between two equations from class (2) is linear with respect 
to u: t = T{t), X = X{t, x), u = U^{t, x)u + U^{t, x), and up to transformations from G^'^^" we 
can assume the coefficient D is not changed. 

Lemma 4. (Ut,Ux) 7^ (0,0) for a local transformation between two equations from class (2) 
only ifDe {ui", e"} mod G^^"'^. 

As an example of discrete equivalence transformations we can give the involution 

t = t, X = —X, u = u-\- ax 
in the couple of equations 

f{x)ut = ie'^Ux)x + ae^Ux and e~°''' f{-x)ut = ie'^Ux)x + ae"-Ux. 
Moreover, this transformation is a discrete invariance transformation for the equation 

g{x)e''^''/^ut = {e''ux)x + ae^n^;, 

iff g is an even function. 

We also investigated some transformations into other classes of reaction-diffusion equations. 
So, using the discrete transformation t = t, x = —x, u = u-\- x/2 we can reduce the equation 

e-^/2^j = e'^iuxx + ul + Ux) 
to the reaction-diffusion equation from the classification of Dorodnitsyn [10]: 

Ut = {e^Ux)x - ^e". 

7 Exact solutions 

We now turn to presentation of some exact solutions for (2). Using our classification with respect 
to all the possible local transformations (i.e. not only with respect to ones from G^'^'^"), at first 
we can obtain solutions of simpler equations (e.g. 6a from Tables 2 or 3) by means of the classical 
Lie-Ovsiannikov algorithm or non-classical methods. Then we transform them to solutions of 
more complicated equations (such as 6b, 6c, . . . ). 

Let us note that the equations with / = 1 are well investigated and the most of exact 
solutions given below were constructed before (see citations in [36]). However, at the best of our 
knowledge, there exist no works containing a systematic study of all the possible Lie reductions 
in this class, as well as the exhaustive consideration of integrability and exact solutions of the 
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Table 6: Reduced ODEs for (27). a ^ 0, s = ±1, 5 = signl 



N 


Subalgebra 


Ansatz u = 




Reduced ODE 


1 




(p{uj) 


t 


V3' = 


2 


(Q4> 


(p{u}) + 21n |x| 


t 


ifi' = 2e'^ 


3 


(Qi) 


ip{ui) 


X 


(ef)" = 


4 


(Q2> 


ip{u}) — In \t\ 


X 


(e^)" = -,5 


5 


(Qi + eQa) 


ip{uj) 


X — et 


(en" = -sv' 


6 


(Q2 +eQ3) 


1/3(0;) — In \t\ 


X — e\n \t\ 


(e^)" = -5(e^' + 1) 


7 


(Qi + eQi) 


ip{oj) + 2et 


— St 

xe 


{e'^y = —suap' + 2e 


8 


{Q2 + aQi) 


ip{uj) + (2a - 1) In \t\ 


x\tr 


{e^y = 5(-aa;¥'' + 2a - 1) 



corresponding reduced equations. That is why we have decided to implement the relevant Lie 
reduction algorithm independently, especially since it is not a difficult problem. 
So, let us consider equation 2.6a 

ut = (e"tx^)^. (27) 

Let us remind that for (27) the basis of A"^^^ is formed by the operators 

Qi = dt, Q2 = tdt - du, = dx, Qi = xd^ + 2du. 

The only non-zero commutators of these operators are [Qi, Q2] = Qi and [Q3, Q4] = Q^. There- 
fore A^'^^ is a realization of the algebra 2A2.1 [34]. All the possible inequivalent (with respect 
to inner automorphisms) one-dimensional subalgebras of 2^2.1 [35] are exhausted by the ones 
listed in Table 6 along with the corresponding ansatzes and the reduced ODEs. 

We succeeded in solving the equations 6.1-6.5. Thus we have the following solutions of (27): 

11,1 1 f~^'^ , ClX + Co 

u = In |cix -I- Col, « = lnl-^-| 

u = (p(x — et) where / d(f = uj + cq. 

J ci-etp 

Using them we can construct solutions for Cases 2.6b-2.6d easily. For example, the transforma- 
tion 2.4 yields the corresponding solutions for the more complicated and interesting equation 

(^e^ + 1)3 ^* = (^""^)^ + ^"^^ 
having the localized density (Case 2.6c): 

ln|ci+c,(e- + 7)|, ^ = l"(- 2,(ei + ^) -f + ^^^)- 
A singular value of the parameter ji for Case 3.6a \s ji = —1. So, the equation 

m (28) 

\ U J X 



u 



ut 

u 



is distinguished by the reduction procedure. It is remarkable that Cases 3.6c and 3.6d are reduced 
exactly to the equation (28). The invariance algebra of (28) is generated by the operators 

Qi = dt, Q2 = tdt + udu, Q3 = dx, Qi = xdx - 2udu 
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Table 7: Reduced ODEs for (28). a / 0, £ = ±1. 



N 


Subalgebra 


Ansatz u = 




Reduced ODE 


1 


(03) 


(p{uj) 


t 


ip' = 


2 


{Qi) 




t 


ifi' = 2 


3 


(Qi> 




X 


(95- V')' = 


4 


(Q2> 


ip{uj)t 


X 


V)' = f 


5 


(Qi + eQ3> 




X — et 


(V" V)' = 


6 


{Q2 + eQs) 


ip{u!)t 


X — eln \t\ 


(93" V)' = -s'p' + 'P 


7 


{Qi + eQ4) 




—et 

xe 


(ip~^ifi')' = —eojifi' — 2s(p 


8 


{Q2 + aQ4) 




x\t\-" 


{ip~^(fi'y = —auxp' + (1 — 2a)ifi 



and is a realization of the algebra 2A2.1 too. The reduced ODEs for are listed in Table 7. After 
integrating Cases 7.1-7.4 we obtain the following solutions of (28): 



2cft 2teoc?e^i^ 

u = coe' , u= 2 — ? — [ ""=7; Z-^P^^ 

cos^ci(x + coj (1 — coe^i-^j^ 

ci £ 2t 

U = -, ;t, U = , U = 



-e + coe'=i(^-^*)' x-et + co' {x + ci)^ + cof^ ' 

Analogously to previous case we obtain by means of transformations 3.5 exact solutions of 
equation 3.6d in the following forms: 



cos2ci(e-^ + co)' (1 -coe'=i<=~")2 ' 



cie"^^ £e-T^ 2te-y''~ 
u= r-rz — 77, u= — — , u 



Another example of equation with localized density is given by Case 3.6f. To look for exact 
solutions for it, at first we reduce it to the equation 3.6a 

ut = {u^U:,)^ . (29) 

As in the previous cases the invariance algebra of (29) 

^'"^^ = {Qi = dt, Q2 = tdt - fi^^udu, Q3 = d^, Q4 = xd^ + 2ii'^udu) 

is a realization of the algebra 2^2.1- The result of reduction (29) under inequivalent subalgebras 
of A^'^ is written down in Table 8. 

For some of the reduced equations we can construct the general solutions. For other ones we 
succeeded to find only particular solutions. These solutions are following: 

1 

, ._L- , , .1 / a (x + co)^ . , if_ 

= Icix + co| ''+1 , u = [cq — £H\X — Et))l^ ^ u=\ V C\\t\ '^+2 



H + 2 2t 
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Table 8: Reduced ODEs for (29). fi ^ 0, -1, a ^ 0, s = ±1, 5 = signt. 



N 


Subalgebra 


Ansatz u = 




Reduced ODE 


1 


(Q3> 




t 


Vs' = 


2 


{Q4) 




t 


<p' = 2At-2(2 + M)<p''+' 


3 


{Qi) 




X 


(¥>>')' = 


4 


{Q2) 




X 




5 


(Qi+eQa) 




X — et 




6 


{Q2 + eQs) 




X — e\n \t\ 




7 


{Qi + eQa) 




xe-'* 


{ip'^ip'y = -sojip' + 2ij,~^£(p 


8 


{Q2 + 0154) 




x\t\-" 





All the results of Tables 7, 8 as well as constructed solutions can be extended to equa- 
tions 3.6b-3.6g using the local equivalence transformations. So for the equation 

ut = (ui'u^)^ + u^u^ (30) 

(e~^ + 7) i+A* 

(Case 3.6f) the transformations 3.7 yield exact solutions in the form 

u = |co(e"^ + 7) - ci|^, u = (cq-I ^ h e'^nt] " |e~^ + 7!"^, 

V e ^ + 7 J 



u=i -[Co ^ I +ci\t\ ^1 |e "" + 71 M+i 

A number of exact solutions were constructed for equations from class (19) (/ = 1) by means 
of nonclassical methods. Starting from them and using local transformations of conditional 
equivalence, we can obtain non-Lie exact solutions for more complicate equations (Cases 6b, 
6c, ...). 

So, N.K. Amcrov [38] and J.R. King [39] suggested to look for solutions of the equation 
Ut = {u^^/'^Ux)x (3.6a, n = —1/2) in the form u = {(p^{x)t + (^^{x))^ where the functions <p^{x) 
and (y5°(x) satisfy the system of ODEs ip].^ = ('^^)^, ^^^x — V^^'P^- ^ particular solution of this 
system is 

16 Ci C2 

x-^ 

And the corresponding solution for the equation 3.6f with /x = —1/2 can be written down as 

u=i6t + c[ + C2e-^)2(e-^ + 7)^. 

8 Conclusion 

In this paper group classification in the class of equations (2) is performed completely. The main 
results on classification are collected in Tables 1-3 where we list inequivalent cases of extensions 
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with the corresponding Lie invariance algebras. Among the presented equations there exist ones 
which have the density / locahzed in the space of x and are invariant with respect to more 
abundant Lie algebras than A^'^'^ . After the tables we write down all the additional equivalence 
transformations reducing some equations from our classification to other ones of simpler forms. 
(In fact, the equations of (2) have been classified with respect to two different equivalence 
relations generated by either the equivalence group or the set of all possible transformations.) 
For a number of equations from the list of the reduced ones we construct optimal systems of 
inequivalent subalgebras, corresponding Lie ansatzes and exact invariant solutions. By means 
of additional equivalence transformations the obtained solutions are transformed to the ones for 
the more interesting and complicated equations with localized densities. 

We describe equivalence transformations within a class of PDEs (2) using the infinitesimal 
and direct methods. The direct method enables us to find the group of all the possible local 
equivalence transformations (i.e. not only continuous ones) in the whole class (2) as well as all 
the conditional equivalence transformations. Moreover, we begin to solve the general equiva- 
lence problem for any pair of equations from class (2) with respect to the local transformations 
(Lemmas 1-4). 

We plan to continue investigations of this subject. For the class under consideration we plan to 
perform classification of potential and nonclassical (conditional) symmetries and finish studying 
all the possible partial equivalence transformations. We also intend to investigate existence, 
localization and asymptotic properties of solutions of initial and boundary- value problems for 
nonlinear convection-diffusion equations. 
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